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Within the QCD factorization approach, we calculate the branching fractions and CP asymme¬ 
try parameters of 12 B —>■ ifo (1430)/('^*^ decay modes under the assumption that the scalar meson 
il'o(1430) is the first excited state or the lowest lying ground state in the quark model. We find 
that the decay modes with the scalar meson emitted, have large branching fractions due to the 
enhancement of large chiral factor r^° . The branching fractions of decays with the vector meson 
emitted, become much smaller owing to the smaller factor r^ . Moreover, the annihilation type 
diagram will induce large uncertainties because of the extra free parameters dealing with the end¬ 
point singularity. For the pure annihilation type decays, our predictions are smaller than that from 
PQCD approach by 2-3 orders of magnitudes. These results will be tested by the ongoing LHCb 
experiment, forthcoming Belle-II experiment and the proposing circular electron-positron collider. 


I. INTRODUCTION 


Although the quark model has made great success in 
describing most of the hadronic states, the lowest lying 
scalar mesons are too light to fit in the quark model. 
Two possible scenarios have been proposed on the ba¬ 
sis of whether the scalars lower than 1 GeV belong to 
the four-quark states or the classical two quark states. 
They are controversial for decades [1]. In scenario-1 

(51) , the scalars such as k( 800), ao(980) and /o(980) are 
naively seen as the lowest lying qq states, and Aro(1430), 
ao(1450), and /o(1500) are the first excited states, cor¬ 
respondingly. In contrast, Arg(1430), ao(1450), /o(1500) 
are treated as the the qq ground states and their first 
excited states are about (2.0 ^ 2.3) GeV in scenario-2 

(52) , in which the lightest scalar mesons are the four- 
quark bound states. 

In the past decade, many B decay modes involving 
scalar mesons have been reported by BaBar, Belle and 
large hadron collider-b (LHCb) experiments. This may 
provide a unique feature to distinguish these two scenar¬ 
ios by the B meson tag [2]. It is thus hoped that the 
combination of the precise experimental measurements 
and the accurate theoretical predictions might provide 
us valuable information on the nature of scalar mesons. 

To achieve this goal, some hadronic Bq{q = u, d, s) de¬ 
cays to scalar mesons have been studied in detail in the 
framework of the QCD factorization (QCDF) [1, 3, 4] 
and the perturbative QCD approach (PQCD) [5, 6]. Us¬ 
ing the QCDF approach, H-Y Cheng et.aZhad calculated 
the branching fractions and direct CP asymmetries of 
most decay modes in refs. [3, 4], such as B -T foK, 
B —iFg (1430)(/)(p) and B —>• ArQ(I430)7r, where most 
results accommodate the data with large uncertainties. 
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Very recently, the decays B —>■ ArQ(1430)Ar(*) dominated 
by 5 —> d penguin operators, have been calculated within 
the PQCD approach [6], however they have not been 
touched in QCDF frame. To complete, we therefore shall 
calculate the branching fractions of H —>■ iFQ(I430)A'^*^ 
decays in QCDF, as well as their CP asymmetry param¬ 
eters. In the experimental side, these observables are 
too small to be measured at current experiments, but 
they are hopeful to be measured in the future experi¬ 
ments such as the updated LHCb experiment, the high 
luminosity Belle-H experiment and the proposing high 
energy circular electron-positron collider. 

The paper is organized as follows. In Section 2, we 
give the analytic formula including the effective Hamil¬ 
tonian, the form factor and all corrections to the ampli¬ 
tudes. Presentation of results and discussions are given 
in Section 3. At last, we summarize this work in Section 
4. 


II. ANALYTIC FORMULA 


In this section, we shall start from the weak effective 
Hamiltonian responsible for B —i 1^0(1430)^1*^*) decays. 
In the standard model, it could be written as [7] 


"Heff = 


Gp 

V2 


VubV:, (C^O'f + C20^) 


10 


~ ( ^2 CiOl -I- C-J^Og^ P CsgOsg) 


i—3 


-l-h.c. (1) 


The explicit form of the operators Oi and the corre¬ 
sponding Wilson Coefhcients Ci at different scale fi could 
be found in ref. [7]. Vu(t)b and are the Cabibbo- 

Kabayashi-Maskawa (CKM) matrix elements. Note that 
Oi _2 are tree operators and others 03 - 10 , 77.89 are pen¬ 
guin ones. 
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In dealing with the nonleptonic charmless B decays, 
the decay amplitude is usually separated into the emis¬ 
sion part and the annihilation part in terms of the struc¬ 
ture of the topological diagrams. According to the fac¬ 
torization approximation based on the heavy quark limit, 
the former part could be written as the product of decay 
constant and form factor. While for the latter one, it is 
always regarded as being power suppressed. In QCDF [ 8 ] 
based on collinear factorization, the contribution of the 
non-perturbative sector is dominated by the form factors, 
and the non-factorizable effect in the hadronic matrix el¬ 
ements is controlled by hard gluon exchange. Thus, the 
hadronic matrix elements of the decay can be written as 

= [ dxTl^{x)<^MAx) 

j io 

+ f f dx [ dyTl\^,x,y)<^B{C)‘^MAx)‘^MAy)^ 

Jo 

( 2 ) 

where Tf- and denote the perturbative short- 

distance interactions, which can be calculated perturba- 
tively. ^x{x) {X = B, Mi, M 2 ) are the universal non- 
perturbative light-cone distribution amplitudes that can 
be estimated by the light-cone QCD sum rules. The form 
factors of B meson to the pseudoscalar (P), the vector 
(y) and scalar [S) mesons transitions being a 

part of the nonperturbative sector of QCD, lack a pre¬ 
cise solution. For the i? —>■ P and P —>■ F transition 
form factors, we will employ the results of the QCD sum 
rule method [9] , since they have been used in calculating 
B —>■ PP, PV and VV modes widely. As the form fac¬ 
tor of as concerned, to the best of our knowledge, 

a number of approaches had been advocated to calcu¬ 
late them, such as QCD sum rule [10], light-cone QCD 
sum rule [11], PQCD [12] and covariant light front quark 
model (cLFQM) [13]. In this work, we use the results 
obtained by cLFQM [13] for keeping consistent with the 
results of other B —>• SP{V) decay modes [3]. 

Following the standard procedure of QCDF approach, 
the emission part of decay amplitude could be written as 


As {BO Ml M2) 

= % E Y.^pbV;^<{y){MiM2mB)F. (3) 

V p—u,c i 


(MiM 2 lOilP )F is the factorizable matrix element, which 
can be factorized into a form factor times a decay con¬ 
stant, as stated before. The effective parameters af can 
be calculated perturbatively, whose expressions are given 

by 


(M 1 M 2 ) 




+ 

+ 


Ci±i CpOis 
Nc 47r 

p!{M2), 


47r^ 


(4) 


with i = 1, ■ ■ ■ , 10. The upper (lower) signs apply when 
i is odd (even), Cp = {N^ — 1)/(2W) with Nc = 3. 
The quantities Vi{M 2 ) account for vertex corrections, 
Pi (Ml M 2 ) for hard spectator interactions with a hard 
gluon exchange between the emitted meson and the spec¬ 
tator quark of the B meson and Pi {M 2 ) for penguin con¬ 
tractions. Similarly, the annihilation contributions are 
described by the terms bi, and 6 i,EWj which have the 
expressions 


(5) 

where the subscripts 1,2,3 of stand for the an¬ 
nihilation amplitudes induced from {V — A){V — A), 
{V — A) {V + A) and {S — P){S + P) operators, respec¬ 
tively. The superscripts i and / refer to gluon emission 
from the initial and final-state quarks, respectively. It 
should be stressed that the decays B and 

B ^ K^*'>+Kq- are only induced by the annihilations 
type diagrams. Hereafter, the order of the a\{MiM 2 ) 
coefficients is dictated by the subscript Mi M 2 , where 
M 2 is the emitted meson and Mi shares the same spec¬ 
tator quark with the B meson. For the bi{MiM 2 ) of 
the annihilation. Mi means the one containing an anti¬ 
quark from the weak vertex, while M 2 contains a quark 
from the weak vertex. The explicit expressions of F, 
Hi,Pi, bi, fcqEW and their inner functions could be found 
in refs.[l, 3]. 

When dealing with the hard-scattering spectator and 
the weak annihilation contributions, we suffer from in¬ 
frared endpoint singularities X = dx/{I — x), that 
cannot be calculated from the first principle in the QCDF 
approach and only be estimated phenomenologically with 
a large uncertainty. Following the arguments of ref. [8], 
we also parameterize these kinds of contributions by the 
complex quantities, Xh and Xa namely, 

Xh = (l + pjie*^-)ln^, (6) 

Xa = {l + PAA^-)hv'^ , (7) 

where A = 0.5GeV. pn^ pa are real parameters, and (pn 
and (pA are free strong phases in the range [—180°, 180°]. 
All the above four parameters should be fixed by the 
experimental data, such as branching fractions and CP 


bi 

b2 

bs 

bi 

^3,EW 

fe4,EW 


Cf ^ .i 

.i 

Cf 
iV2 L 
Cf 
iV2 L 
Cf 
iV2 L 
Cf 
N' i 


C 3 AI + C5(A* +aI) + NeCeAl 

CiA\ + CeA(\ , 

+ < 77(^3 + A^) + NcCs,A\ 
[CioAl + CgA* ] , 
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asymmetries. In the so-called PQCD approach [14], one Including the emission and annihilation contributions, 
can eliminate these end-point singularities by keeping all the decay amplitude can be finally given as 
small transverse momenta of gluons and inner quarks. 


A{B- 


K-K*°) 


^ p—u,c ^ '' 


+ /B/iC*/ic(^2^S + ^3 + &3.Ew)j 


( 8 ) 


Gf 


AiB-^K*-K^) = ^ E - 


p—u,c 


P P 

“4 - 2^10 


- « - ) fKFo^^° - m^.) 

^ / K*-K0 


+ /B/if5/if + ^3 + ^3,Ew);x* 


K^-KO f’ 


(9) 


A{B° Ko°) = ^ E («4 - ^a?o - K - - mj^) 

^ p—u,c \ ^ y K Kq 

+ fBfKS /k [ (^bs + b4- - (foa^EW + fe4,Ew)^ 


ITK^O 


+ 64 - -b. 


4 — :re4,EW 


( 10 ) 


a{W^Xk°) = 


Gf 

V2 


E I “ (“4 ” ^“1° “ (.rnK){ml - m\,) 

}—U,C K ^ ^ ^0 

+ fBfK*fK [ (^bs + b4- -(63 ,EW + & 4 .Ew)^ 


_ -I- I &4 — ■;^&4,EW _ 

K^Oro \ Z y KOK' 


( 11 ) 


A{B°^Kq~''K ) = E + ^4 + &4.Ew)e'* + E-- + - 2^4,Ew'j (12) 

A{B^^K^Kl~) = ^ AG)|/5/^./^[(6;^(5p-I- 64 -I-64_ew)e'+;x*-+ f&4 - 2^4,Ew'j _ ]p (18) 

^ P—U,C ^ ) 


A{B- ^ K*-K*^) = ^ ^ aG)| - U - iafo + (a| - 

’ p—u,c \ ^ 


X 


K-K*o 


2/if*^o^ (m^-)wBPc 


— /bIk* /if* (^2^S + ^3 + fe3.Ew) 


if-if*o f’ 


(14) 


il(B K*°) = % E («4 - ^<0 - (^6 - ^aD^’x *) ‘2fK-F^^^°{mj^.)mBPc 

^ P—U,C I, ^ ^0 

-fBfK*fK* (&2^S + ^3 + ^3,Ew);x*-if*o f ’ 


(15) 
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Adf^ir^AT) = % E - 


•/2 


p—u,c 


4-o«?0 + K-oW“ 




-/s/if^/if* [ ( ^3 + ^4 - ■;^(^3,EW + &4,Ew) ) + ( ^4 - :^&4,EW ) ], 

V ^ \ ^ J K*°K'° 


(16) 


~*0 'y '^ ^(d) J f _ ^„P _ („P _ ^„P\,r,P 


A{W^KoK*°) = 


V2 


p—u,c 


< - ^<0 - i< - ) 2/^.Ff (m?,.)mBPc 


/ Kg°K*o 


-/s/ETj/if* [ (^bs + 64 - -(fea.EW + & 4 ,Ew)^ 


_ + I &4 — T^fe4,EW _ 

K*°K*o V ^ / ifOK* 


1 , (17) 


a(b° a k;*k>-) = ^ E 

^ p—U 


, “ /B//f*/Er[ (^l^S + ^4 + &4.Ew)if*+if.- + ^^4 - 2^4,EW^ 


>,(18) 


if—_fc: 


A(i?° ^ K*+K*-) = ^ 


E -/B/if-/if[(MS + &4 + 64 ,Ew)x* 

p—u,c \ 


+ if*- + ^4 - ■r&4, 


EW 


>,(19) 


K*-K- + 


r 


where 


rf(/4) = 

'■f (m) = 


2m\ 


mb{^l)[mu{^J) + TOs(/i)] ’ 


2m 


ifn* 


„if* 


(m) = 


mb(/i)[mq(/i) - ms(/i)]' 
2 tok. /j^. (^) 
rribifi) Ik* 


( 20 ) 


of Xq( 1430) have the signs flipped from SI to S2. The 
definition of the form factors for the B ^ S transitions 
are given by [13] 


{S{p')\A,\B{p)) = -^ 


Pp- 


nriQ — Trig 


^2 ^0 W ) 


Qp FrHr) 


(23) 


III. NUMERICAL RESULTS AND DISCUSSION 


To proceed, we shall present numerical results obtained 
from the formulas given in the previous section. Firstly, 
we introduce the adopted parameters in our calculation. 
Secondly, we give the numerical results and show the the¬ 
oretical errors due to uncertainty of some parameters. At 
last, some discussions and comparisons will be added. 

The scalar meson, unlike the pseudoscalar one, has two 
kinds of decay constants, namely the vector decay con¬ 
stant fs and the scale-dependent scalar decay constant 
/s, the definitions of which are give by: 


{Koip)\s'lpd\0) = fK-Pp., 

{Kl{p)\sd\0) = mK-jK-. 

The two decay constants satisfy 

ms{p) - md(,p) j 

JK* = - JK*, 

° rriKi ° 


( 21 ) 


( 22 ) 


where nig^p) and rridip) are the running current quark 
masses. It should be stressed that the decay constants 


where = {p + p')^ and = [p - p')^. In cLFQM , 
the momentum dependence of the form factor could be 
parameterized in a di-pole model form [13], 

PC 2,, ^ _ Pj^) _ C24i 

> l-a{qyml) + biqym%y ^ > 

Together with the decay constants, the parameters F{0), 
a and b for B ^ S transitions in the different scenarios 
are summarized in Table 1. 

The twist-2 light-cone distribution amplitude (LCDA) 
$s(a:) and twist-3 4>g(a;) and $g(a;) for the scalar meson 
Kq made of the quarks sd are given by 


{K;{p)\s{z2hpd{zi)\0) = pp r dxed-P-^^+^P-^^'><^>s{x), 

Jo 

(K*o{pMz2)d{zO\0) = [ dxed-P---+-P-^^)^yx), 

{Koip)\s{z2)<Jp„d{zi)\Q) 

=-mj^^ippz, - P.zp) (25) 
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TABLE I: The Parameters of the Js'o(1430) in different scenarios 



(0) 

a 

b 


a 

b 

/^.(MeV) 

Bi 

Ba 

SI 

0.26 

0.44 

0.05 

0.26 

1.52 

0.64 

-300 ± 30 

0.39 ±0.05 

-0.70 ± 0.05 

S2 

0.21 

0.44 

0.05 

0.21 

1.52 

0.64 

445 ± 50 

-0.39 ± 0.09 

-0.25 ±0.13 


TABLE II: Summary of input parameters 



A A 

Ru 7 rgo 

Xb- 

As Qe 




0.225 0.818 

0.376 68° 250MeV 1.54ps 

1.67ps 

0.35 1/132 



/s 

niB fK 

fK* fk* m-K 

rriK* 

niK* 

pB^K 

aB^k^ 

^0 

236MeV 

5.36GeV 131MeV 

221MeV 175MeV 0.49GeV 

0.89GeV 1.43GeV 

0.35 

0.34 


with z = Z 2 — zi, X = 1 — and their normalizations are 


[ dx<i>s{x) = fs, 
Jo 


dx^six) = / dx^g{x) = fs- 


(26) 

(27) 


Jo Jo 

The above definitions of LCDAs can be combined into a 
single matrix element as 


{K;{p)\s2p{z2)dUz,)\0) = i 

4 Jo 

X |/i$s(a:) + ms “ (x^i,p^z‘' ^ | 

“(28) 

In general, the twist-2 light-cone distribution amplitude 
<i)s has the form 


$s(x, //) = /5 (m) 6x(1 - x) 


X 


Bo{^i)+Y.BMC^\2x-l) 

m—1 


(29) 


where Bm (/r) are scale-dependence Gegenbauer moments 
and Cm'^{u) are the Gegenbauer polynomials. The Bm 
of different scenarios are also presented in Table I. As for 
the twist-3 distribution amplitudes, we shall adopt the 
asymptotic form for simplicity, shown as 

^s{x) = fs, $5(x) =/s6x(l - x). (30) 

For the LGDAs of K^*\ we will employ the formulaes 
obtained from the QCD sum rules [15, 16]. The other 


used parameters, such as the CKM elements, the decay 
constants of the pseudo-scalar and the vector, and the 
form factors of B —>■ K^*\ are also listed in Table II for 
convenience. 

Now, we turn to discuss the numerical results of the 
concerned decay modes. As we had stated in previous 
section, when calculating the hard spectator and the an¬ 
nihilation contributions, two endpoint singularities are 
parameterized phenomenologically by four free parame¬ 
ters, namely ph,4>h and pA,PA, which should be deter¬ 
mined from the experimental data. In ref. [3] , a global fit 
of PA and (/)A to the B —>• SB data indicates pA = 0.15 
and 4>a = 82° with = 8.3, so in this work for the 
central values (or “default” results), pa,h = 0.15 and 
4>a,h = 82° are adopted. In Table III, we list the cal¬ 
culated branching fractions of i? —>■ ArQ(1430)iT^*^ de¬ 
cays, where the first uncertainty is due to the variations 
of i?i _3 and /s, the second comes from the form factors 
and the strange quark mass, and the last one is induced 
by weak annihilation and hard spectator interactions in 
ranges pa,h G [0, 0.3] and (I>a,h & [0,180°]. 

Here, we shall take B~ -A 7F“7Fg°(1430) and B~ -A 
iFg” (1430)iF° as examples to illustrate the relative size of 
each contribution. The decay amplitude formula of each 
mode has been presented in eqs.(8) and (9), respectively. 
The first decay mode is characterized by i? —>■ P transi¬ 
tion with scalar meson emitted; while the second decay 
mode is characterized by P —^ S' transition with pseu¬ 
doscalar meson emitted. Because of the small vector de¬ 
cay constant of scalar meson A'g(1430), 

is suppressed relative to (ixix)- Therefore, the 

decay width of the first channel should be suppressed 
comparing with the second one. The numerical results of 
these two decay modes are given as below: 


A{B- -A K-K*°) oc VubV*a (0.53 -f 0.14z) -h VcbV^^ (0.58 -h 0.081) -f (-0-03 -h O.OH) -h VcbV;^ (-0.03 -h O.OH), 


emission diagrams 


annihilation diagrams 


(31) 
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TABLE III: The branching ratios of B ^ Aq (1430)A'(*^ in the unite 10 


Decay Mode 


SI 


S2 


b" AS“(1430)A:° 

B- Ko-(U30)K° 
B° ^K°KS°{U30) 
B- a:-a:5°(1430) 
B° A'*+(1430)A:- 
# ^ A+A*~(1430) 
B*" A;“(1430)A'*° 

B” Ao-(1430)A*° 
B° A*V5°(1430) 
B- A*-Ao°(1430) 
'W’ A*+(1430)A*- 

-» A*+A*~(1430) 


I I -19+3.80+7.07+4.91 
99-3.57-5.05 
r ofj + l.73+5.13+1.73 
^•OU-i 37-2.36-2.30 
Q9 Q'T+0.41+7.80+5.49 

OZ.Z/ _7 g3_5 4g_g 30 

9Q 'T'l+6.67+6.73+2.61 
AO. /1-5.60-4.61-3.64 
rj Q7+0.43+0.01+0.63 
U-yZ-o.31-0.01-0.44 
+3.31+0.07+7.28 
-2.55-0.10-4.83 


9 00+1.20+1.95+2.67 
Z.oy_o.85-0.90-2.00 
n -,.+0.54+1.40+1.17 
1-1^-0.38-0.56-0.92 
An /17+13.36+6.09+6.06 
_10.77-5.38-6.16 
QQ 70+10.33+5.52+3.37 
OO. /U_g 47_4 g2_3.94 
n CQ+0.45+0.02+0.14 
U.OO_o.29-0.03-0.05 


8.33 


1.07 


+0.72+0.03+2.27 

-0.47-0.04-0.97 


U.UO_o.02-0.04-0.06 
Q 02+0.08+0.21+0.05 


^-0.08-0.18-0.04 

+2.54+0.25+0.53 


10.86_2 24-0.24-0.75 

-,9 7 -,+3.14+0.26+1.16 
IZ. Z 1-2.72-0.26-1.07 
r» QQ+0.14+0.00+0.22 

U.OO_o. 10-0.00-0.16 

14 c:.^+5.75+0.00+12.48 
l^-^^-4.44-0.00-g.31 


ri ^iQ+0.03+0.02+0.03 
U.UO_o. 00 - 0 . 00 - 0.02 
r. -, ^+0.06+0.11+0.03 

1^-0.05-0.08-0.03 
on 1 q+5.35+0.50+0.42 
ZU.10_4 49 _ 0 . 49 - 0.59 
91 71 + 5 . 51 + 0 . 54 + 0.30 
Zl. Z 1-4.66-0.53-0.24 
r. 1 1 +0.09+0.00+0.09 

U.11-0.06-0.00-0.06 

o/, + l.26+0.00+3.94 
^•Q^-0.83-0.00-1.75 


A{B- K*-K°) oc VubV*a (0.20 + 0.05t) + VebV;*^ (0.22 + 0.03t) + (0-08 - O.OH) + VebKd (-0-04 - O.Oli). 

V V V 

^ V V 

emission diagram annihilation diagrams 

(32) 


From these two equations, it is apparent that the emis¬ 
sion diagrams are dominant. However, there is a large 
enhancement from Oe.s operators due to the fact that the 
chiral factor r^° = 12.3 at /r = 4.2 GeV is much larger 
than = 1.5 owing to the larger mass of iCo(1430). 

It follows that (og — ^a.g)rx° is much greater than 
(ttg — and (ug — • This compensate with 

the suppression of scalar meson decay constant to result 
in a larger branching ratio of B~ —>■ A"“iCg°(1430) with 
scalar meson emitted. 

In the B —>• PPiy) and VV decay modes, the weak 
annihilation contribution is usually expected to be very 
small because it belongs to the next leading power correc¬ 
tion. However, one can see from Table HI that the uncer¬ 
tainties induced by the weak annihilation are very large, 
even much larger than the central values in some decay 
modes, such as B —>■ Kq (1430)A"° and pure annihila¬ 
tion mode B —?> A"Q^(I430)Ar*“. This phenomenon can 
be understood as follow, when discussing the penguin- 
induced annihilation diagram of B —^ PP mode, the de¬ 
cay amplitude is helicity suppressed heavily because the 
helicity of one of the final-state mesons cannot match 
with that of its own quarks. On the contrary, this kind 
of helicity suppression can be alleviated when the scalar 
meson involved due to the nonzero orbital angular mo¬ 
mentum Lz of the scalar meson. 

From Table HI, one could notice that our predicted 
central values for the branching ratios of Kq (I430)A'° 
and Arg“(I430)A''^ based on S2 are smaller than the re¬ 
sults based on SI by a factor of 4 ^ 5. In contrast, 
the central values of of Ar°A'g°(1430) and K~Kq^[ lAAQ) 


based on SI are a bit smaller than those of S2. For 
Ar*°Arg°(1430) and AT*”Arg°(1430), the predicted cen¬ 
tral values of S2 is larger than those of SI by a factor 
2. Although there is difference between the central val¬ 
ues of different scenarios, we cannot distinguish two sce¬ 
narios due to large uncertainties taken by annihilation 
diagrams, unless one approach were proposed to deal 
with annihilations effectively in QCDF. In previous stud¬ 
ies [3], by comparing the calculated branching ratios of 
B —^ KqTt and B K^cj) with experimental data, one 
concluded that the S2 is prefered, i.e. Kq is very likely 
the lowest lying state. If so, the branching fractions 
of Ar'^A'g°(I430) and A"“Arg°(I430) could be measured 
by analyzing the data of AT AT+tt” and K~K~^Tr~. Un¬ 
fortunately, the predicted results of ATg (1430)Ar*° and 
ArQ“(1430)Ar*° is too small to be measured with the cur¬ 
rent data of Belle, but they are hopeful to be measured 
in the forthcoming Belle-H by analyzing the four-body 
final states of AT+AT^Tr+Tr”. 

Within the framework of PQCD, X. Liu et.al had cal¬ 
culated these decays [6], and both their results and ours 
are below the experimental upper limits. Comparing our 
predictions and theirs, we find that for pure annihilations 
they obtained rather larger branching fractions by keep¬ 
ing the transverse momenta, which are larger than ours 
by more than two order of magnitude. For decays with 
emission diagrams, they also got the larger branching ra¬ 
tios with rather large nonfactorizaion diagrams based on 
both two different scenarios. We hope the future mea¬ 
surements could distinguish two different frameworks. 

The CP asymmetries have not been observed in any B 
decays involving a scalar meson. For the charged decay 
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TABLE IV: The direct CP asymmetries (%) of B A:5(1430)A'(*) 


Decay Mode 

SI 

S2 

B- A*-(1430)A:° 
B- A-A'o°(1430) 

r 97+0.50+0.31+2.59 
i -0.59-0.57-2.82 

9 nA+0-^0+0-85+5.46 
^•119-0.65-0.69-6.67 

99 tti+4.90+5.63+19.61 
ZZ.01_7 57_g 3g_22.86 

9 f?n+l-61+0-59+3.52 
Z.UU-i 76-0.59-5.47 

B- K*-{U30)K*° 
B- A'*-A'o°(1430) 

1ft on+0.94+0.65+1.27 
10.0U_;L. 12-0.89-0.99 
r j^o + l.30+0.17+11.40 
9.UO i3y o. 17-13.52 

qi 09+4-67+4.28+5.06 

Ol.UZ_g 5g_7 g5_3 go 

O 0/1+2-54+0.23+7.57 
U.O^2.64-0.22-9.22 


modes, according the definition, 


Sf appearing in above equations are defined by 


jdir _ MB f) - A{B+ -» /) 

CP a{b-^ f) + A{B+^ f-y ^ > 

the predictions of the direct CP asymmetries based on 
two different scenarios of scalar mesons are summarized 
in Table IV. The resource of this asymmetry is the inter¬ 
ference between the emission diagrams and annihilations. 
The large uncertainties in the latter lead to the large er¬ 
rors in the direct CP asymmetries, as shown in the table. 

For the neutral decay modes, the situation becomes 
more complicated due to the fact that both B^ and B 
could decay to the same final states. For illustration, we 
will take B^{B^) —>■ /TQ°(1430)iC° as an examples. Four 
decay amplitudes, Aj, Aj, Af and Ajr are denoted as 


Af = {K*olM\By, Af = {KlKyBM 

A/ = {K*K°\^), Af = (k;k°\B°). 

Then, we can write down the CP asymmetry as 

\AM + \AM-\Af\^-\Af\^ 

\Af\^ + \Af\^ + \Af\^ + \Af\^- 


(34) 


(35) 


Due to B° — B° mixing, the four time-dependent decay 
widths are given by (/ = /Tg°(1430)/T°) 


T{B°{t) ^ /) 

r(B°(<) ^ /) 

T{B°{t) ^ /) 

r(B°(<) ^ /) 


[1 + Cf cos Amt — 5/ sin Amt], 

e-^^l{\AfM + \AfM 

[1 — Cj cos Amt + Sf sin Amt], 

e-^^^{\AfM + \AfM 

[1 + Cj cos Amt — Sf sin Amt], 

e-^^l{\AM + \AM) 

[1 — Cf cos Amt + Sf sin Amt] ,(36) 


where Am stands for the mass difference of two mass 
eigenstate of B°/B meson, and F for the average decay 
width of the B meson. The auxiliary parameters C / and 


Cf = 

\AM-\AM 

I+P + I+P’ 

(37) 

Sf = 

2Im(A/) 

(38) 

1 -b \Af/Af\'^' 

^f = 

Vtbv;, M 
yivtd Af^ 

(39) 


Replacing / by /, we could obtain the formulaes for Cf 
and Sf^ correspondingly. If the experiment could find 
the values of Cy, Sf, Cf and Sf, we thus can obtain four 
new parameters: 

C = \{Cf + Cf), AC = \{Cf-Cf), (40) 

S = \{Sf + Sf), AS=\{Sf-Sf). (41) 

Physically, S is the mixing-induced CP asymmetry and 
C is from the direct CP asymmetry, while AC and AS' 
are CP-even under CP transformation A/ —>■ 1/Aj. In 
Table. V, we present our estimations of Aqp, C, AC, S 
and AS for the final states Kq K^, Kq K*^, Kq'^K~ 
and KMK*~, under two different scenarios. It should 
be stressed that in PQCD [6], there is no CP asymme¬ 
tries because of the absence of tree operators. However, 
in QCDF, by including the penguin contractions (Pi), 
the charming penguin namely, an extra strong phase to¬ 
gether with another one from annihilations might lead to 
the large CP asymmetry, as shown in the table. Again, 
the large uncertainties of annihilations result in the large 
error for the pure annihilation decay modes, especially 
under S2. We hope these parameters can be measured in 
future colliders, such as high luminosity Belle-II, LHC-b 
and even higher energy e“'"e“ collider. 


IV. SUMMARY 

In this work, we studied the B —>■ Kq{ 14 S 0 )K^A decays 
under two different scalar meson scenarios by using the 
QCD factorization approach. We calculated the branch¬ 
ing fractions and the CP asymmetry parameters. It is 
found that the decay modes with the scalar meson emit¬ 
ted, have large branching fractions due to the enhance¬ 
ment of large chiral factor r^°, with some of the branch¬ 
ing fractions around the corner of Belle II. In contrast. 















TABLE V: The CP-violating parameters Acp, C, AC, S and AS (%) of B KS{U30)K^*\ 


Decay Mode 

Scenario 

Acp 

c 

AC 

s 

AS 


SI 

c, r^+0.04+0.12+0.08 
U.OZ o.o4-0.10-0.13 

n nc^+O-OO+o.oo+o.oi 
U.UO o. 00 - 0 . 00 - 0.01 

^ p>q+ 0 . 00 + 0 . 01 + 0.00 

U.UO o. 00 - 0 . 01 - 0.00 

I pjpj+ 0 . 00 + 0 . 01 + 0.00 

i.UU- 0 . 00 - 0 . 01 - 0.00 

n AA+o. 00 + 0 . 00 + 0.00 
U.UU-o. 00 - 0 . 00 - 0.00 

S2 

0 Qri”l"0'03+0.06+0.08 
—u.yu 0.02-0.04-0.08 

ri ^r+0.01+0.01+0.03 

U.UO o. 00 - 0 . 01 - 0.01 

r» p\ 9 + 0 . 01 + 0 . 01 + 0.02 

— U.UZ Q QQ Q QJ 0 QQ 

1 ^p\+ 0 . 00 + 0 . 00 + 0.00 

i-UU- 0 . 00 - 0 . 00 - 0.00 

PI p. 9 + 0 . 00 + 0 . 00 + 0.00 

U.UZ- 0 , 00 - 0 . 00 - 0.00 


SI 

^ Qo+0.02+0.01+0.02 

—u.yo o. 00 - 0 . 01 - 0.01 

pj qq+0.01+0.00+0.14 
t».OO-0.07-0.03-0.08 

p. 9PJ+0.01+0.00+0.14 
^•^^-0.06-0.03-0.08 

ri Q-1+0.03+0.02+0.04 
U.yi-o.03-0.02-0.22 

p\ P 19 + 0 .03+0.03+0.03 

U.UZ o 02 - 0 . 02 - 0.22 

S2 

I pjpj+ 0 . 00 + 0 . 00 + 0.00 

-L.UU-o. 00 - 0 . 00 - 0.00 

ri 9’^+0.09+0.08+0.21 
-0.24-0.23-0.33 

ri "I o"l"0.09+0.08+0.20 
'^•-*-®-0.24-0.23-0.33 

n Of-+0.07+0.06+0.07 

U.OD o. 01 - 0 . 00 - 0.22 

n 1 q+0.07+0.06+0.07 
U.-LO-o. 01 - 0 . 00 - 0.22 

K*+K- 

SI 

r» '^Q+ 0 . 01 + 0 . 00 + 0.02 

'^-0.01-0.00-0.05 

r, rj^+0.01+0.01+0.04 
01-0.00-0.04 

U.UO o.02-0.00-0.06 

r» + 0 . 02 + 0 . 01 + 0.02 

'-'•^^-0.02-0.01-0.05 

PI o/j+0.03+0.01+0.10 
'-'•0^-0.03-0.01-0.04 

S2 

rj 9Q+0.06+0.01+0.35 
y-'^y-0.01-0.00-0.98 

r, r|9+0.03+0.01+0.27 

U.UZ Q Q 2 o.oO- 0.10 

rj 9r+0.04+0.01+0.29 
U.ZO o.03-0.01-0.26 

n qc+o.05+0.04+0.50 
U.OO o.05-0.02-0.11 

PI 9 PI+ 0 . 12 + 0 . 01 + 1.02 

— U.ZU o. 10-0.02-0.35 

k;+k*- 

SI 

n Qc+o. 00 + 0 . 00 + 0.00 
u.yo o. 00 - 0 . 00 - 0.01 

^ + 0 . 00 + 0 . 00 + 0.01 
U.U-L- 0 . 00 - 0 . 00 - 0.02 

ri ^9+0.01+0.01+0.03 

U.UZ o. 00 - 0 . 00 - 0.02 

p. PJ-,+ 0 . 00 + 0 . 00 + 0.00 

U.Ul- 0 . 00 - 0 . 00 - 0.00 

PI QQ+ 0 . 00 + 0 . 00 + 0.02 

—u.yo o. 00 - 0 . 00 - 0.01 

S2 

ri oq”I"0'02+0.00+0.05 

U.OO-0.00-0.00-0.66 

pj pj-i +0.01+0.00+0.14 

y.ui-o.01-0.00-0.02 

nil ■l"0.03+0.01+0.23 

y-j-i-o.02-0.00-0.11 

rj PJ9+0.01+0.00+0.48 

U.UZ o 02-0.00-0.02 

PI Q-j+0.02+0.00+0.63 

— U.yi o.02-0.00-0.05 


the branching fractions of decay modes with the vector 
meson emitted, are much smaller. Moreover, the anni¬ 
hilation contributions take large uncertainties because of 
the free parameter from endpoint singularity. For the 
pure annihilation type decays, we predicted very small 
branching fractions, which are 2-3 orders of magnitudes 
smaller than the results from PQCD. Some of the decay 
channels are hopeful to be measured in future colliders, 
such as Bellc-II, LHC-b and even high energy e+e” col¬ 
lider. 
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